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ABSTRACT 

The  rapid  cooling  of  the  mixed  layer  in  the  Arabian  Sea  during  the  southwest 
monsoon  onset  is  examined.  The  atmospheric  forcing  fields  over  the  Arabian  Sea 
during  the  onsets  in  1979,  1984  and  1985  are  analyzed.  Realistic  wind  speed  and 
evaporative  heat  flux  increases  occur  soon  after  onset.  The  solar  radiation  flux  anal- 
yses suggest  that  both  the  Navy  Operational  Global  Atmospheric  Prediction  System 
and  the  National  Aeronautics  and  Space  Administration  four-dimensional  data  assimi- 
lation models  are  too  simplistic  for  accurate  solar  radiation  flux  calculations  in  the 
tropics.  The  solar  radiation  fluxes  appear  to  be  too  low  during  the  pre-onset  periods 
and  too  high  during  the  post-onset  periods.  The  atmospheric  forcing  fields  are  then 
applied  to  a  one-dimensional,  ocean  mixed  layer  model  at  four  locations.  The  magni- 
tude of  the  predicted  sea-surface  temperature  decrease  is  0.5  to  1.0°C  during  the  1984 
and  1985  monsoon  onsets.  The  results  suggest  that  local  atmospheric  forcing  accounts 
for  most  of  the  mixed  layer  cooling  on  a  time  scale  of  less  than  ten  days  and  away 
from  any  major  ocean  current  regions.  The  direct  heat  loss  to  the  atmosphere  contrib- 
utes more  to  cooling  of  the  mixed  layer  than  does  entrainment  mixing  at  the  base  of 
the  layer. 
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I.  INTRODUCTION 

A.  PURPOSE  AND  HYPOTHESIS 

The  purpose  of  this  study  is  to  investigate  the  upper  ocean  response  to  atmos- 
pheric forcing  that  is  associated  with  the  onset  of  the  southwest  monsoon  in  the 
Arabian  Sea.  Estimates  of  atmospheric  surface  fluxes  as  the  monsoon  begins  are  then 
used  with  the  Garwood  (1977)  one-dimensional  bulk  mixed  layer  model  to  simulate  the 
effects  on  sea-surface  temperature  (SST)  and  mixed  layer  depth  as  the  monsoon  begins. 
A  second  objective  of  this  study  is  to  determine  which  atmospheric  forcing  parameters 
can  account  for  the  observed  cooling  of  the  Arabian  Sea  immediately  following  onset 
of  the  summer  monsoon.  The  SST  decreases  rapidly  immediately  after  onset 
(Krishnamurti,  1981),  even  though  the  solar  heat  flux  would  be  expected  to  lead  to 
increasing  temperatures  during  this  season.  The  hypothesis  of  this  study  is  that  the 
increased  wind- generated  mixing  associated  with  the  onset  of  the  monsoon  is  respon- 
sible for  the  immediate  cooling  at  the  ocean  surface.  This  hypothesis  is  only  applicable 
away  from  the  region  of  the  African  coast,  where  upwelling  is  likely  to  be  the  dominant 
cooling  mechanism,  and  on  the  10-day  time  scale  of  monsoon  onset. 

One  motivation  for  this  study  is  to  identify  the  causes  for  the  Fleet  Numerical 
Oceanography  Center  (FNOC)  upper-ocean  model's  predictions  of  excessive  cooling  of 
the  upper  ocean  after  the  southwest  monsoon  is  established  (personal  communication, 
M.  Clancy,  FNOC).  Understanding  the  various  ocean  cooling  mechanisms  may 
contribute  to  improved  predictions  and  analyses  of  upper  ocean  thermal  structures  in 
this  area.  Some  studies  have  also  suggested  a  possible  link  between  SST  anomalies  in 
the  Arabian  Sea  after  onset  and  rainfall  amounts  in  India  (Babu  et  al.  1985).  A  better 
knowledge  of  the  atmospheric  forcing  that  leads  to  the  cooling  of  the  Arabian  Sea  may 
then  lead  to  better  predictions  of  rainfall. 

B.  BACKGROUND 

The  southwest  monsoon  over  the  Arabian  Sea  occurs  annually  during  late  May 
or  early  June.  It  is  characterized  by  a  dramatic  shift  of  surface  winds  from  weak  east- 
erlies or  westerlies  to  strong  southwesterly  winds.  Krishnamurti  et  al.  (1981)  found 
that  a  tropical  cyclone-like  feature,  termed  the  onset  vortex,  accompanies  the  onset  of 
the  monsoon  in  about  60%  of  the  monsoons  observed  during  a  68-year  period.    The 


onset  vortex  normally  forms  off  the  southwest  coast  of  India  at  about  10°N,  70°E  and 
moves  slowly  northward  and  then  westward  to  the  Saudi  Arabian  peninsula  area.  The 
strong  westerlies  to  the  south  of  the  onset  vortex  help  to  establish  the  southwest 
monsoon.  These  monsoonal  winds  bring  rain  to  the  western  coast  of  India  and  areas 
northward  and  inland.  In  addition  to  stronger  winds,  the  onset  vortex  also  reduces  the 
solar  insolation  at  the  surface  due  to  increased  cloud  cover.  Both  of  these  conditions 
tend  to  reduce  the  SST  during  the  period  of  the  monsoon.  Briefly,  the  mixed  layer  heat 
budget  is  represented  by: 

QML  =  Qsfc  +  Qent  +  Qadv  +  Qupw  <M) 

where  Qjyn  is  the  change  in  heat  content  in  the  mixed  layer  ;  Qs^c  is  the  surface  heat 
flux;  Qent  is  the  entrainment  heat  flux  at  the  base  of  the  layer;  Qacjv  is  the  heat  change 
due  to  horizontal  advection;  and  Qumv  is  the  heat  flux  due  to  upwelling.  These  terms 
are  represented  in  Fig.  1.1.   The  heat  budget  is  discussed  in  detail  in  Chapter  II. 

Krishnamurti  (1981)  found  that  the  SST  initially  decreased  by  a  maximum  of 
2-3°C  in  a  5-6  day  period  soon  after  the  wind  strengthened  in  1979.  He  suggested  five 
possible  mechanisms  for  the  cooling  of  the  Arabian  Sea  following  onset: 

(1)  Diminished  solar  radiation  due  to  increased  cloud  cover; 

(2)  Northward  flux  of  heat  due  to  cold  oceanic  eddies  to  the  north; 

(3)  Coastal  upwelling  and  downstream  shedding  of  cold  eddies; 

(4)  Strong  evaporation  in  the  region  of  strong  winds;  and 

(5)  Strong  wind  stress-induced  upwelling. 

This  study  will  also  examine  the  mechanism  of  turbulent  mixing  in  the  upper  ocean 
and  the  associated  entrainment  of  cooler  water  from  below.  Each  of  these  mechanisms 
are  discussed  in  detail  in  Chapter  II. 

C.       DATA  SOURCES 

Two  four-dimensional  data  assimilation  models  are  used  as  the  sources  of  atmos- 
pheric forcing.  This  type  of  data  is  used  in  place  of  direct  observational  analyses  to 
obtain  complete  and  consistent  sets  of  surface  fluxes  on  the  space  and  time  scales 
required  for  ocean  prediction.  The  observation  network  in  the  Arabian  Sea  is  insuffi- 
cient, even  during  the  First  Global  atmospheric  research  program  Global  Experiment 
(FGGE),  to  provide  the  required  data  coverage.  The  data  assimilation  technique  of 
both  data  sources  is  to  use  predicted  fields  to  fill  in  the  gaps  where  no  observations 
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exist.  The  National  Aeronautics  and  Space  Administration  (NASA)  four-dimensional 
data  assimilation  model  is  the  source  of  the  forcing  during  1979.  These  data  were 
obtained  during  the  FGGE  Special  Observation  Period  (SOP)  II.  The  source  for  1984 
and  1985  is  the  Navy  Operational  Global  Atmospheric  Prediction  System  (NOGAPS), 
which  is  maintained  by  FNOC.  This  data  source  is  described  in  more  detail  in 
Appendix  A. 

Chapter  II  describes  the  atmospheric  forcing  fields  from  the  1979,  1984  and  1985 
summer  monsoon  seasons.  A  comparison  is  made  between  the  atmospheric  forcing 
fields  from  each  data  source  and  interannual  differences  are  noted.  The  1979  data 
should  be  more  accurate  since  FGGE  provided  more  direct  observations,  but  no 
concurrent  ocean  thermal  structure  data  are  available  to  verify  the  ocean  mixed  layer 
predictions.  The  1984  and  1985  data  are  based  on  fewer  observations,  but  the  SST 
analyses  used  by  NOGAPS  are  concurrent  with  the  atmospheric  forcing.  Therefore, 
data  for  all  three  years  are  analyzed  to  ascertain  the  accuracy  of  each  data  set  in 
comparison  with  previous  studies  listed  in  Table  1. 

D.      ATMOSPHERIC  FORCING  OF  MIXED  LAYER  MODEL 

The  second  part  of  the  study  involves  the  application  of  the  atmospheric  forcing 
parameters  to  the  Garwood  ocean  mixed  layer  model.  This  model  assumes  no  hori- 
zontal or  vertical  advection  so  that  a  local  heat  balance  is  assumed.  It  is  initialized 
with  climatological  temperature  and  salinity  profiles  at  one  meter  intervals  in  the  upper 
200  m.  The  model  integrates  hourly  values  of  net  surface  heat  flux,  solar  radiative  flux 
and  wind  stress,  and  predicts  the  changes  in  the  mixed  layer  temperature  and  depth. 
Specifically,  this  study  utilizes  the  predicted  maximum  daily  mixed  layer  depth  and 
associated  temperature,  the  minimum  daily  mixed  layer  depth  and  associated  tempera- 
ture, and  the  mean  daily  mixed  layer  depth  and  associated  temperature. 

A  major  limitation  to  this  study  is  the  lack  of  verification  SST  data.  Therefore, 
the  main  emphasis  will  be  to  determine  the  relative  effects  of  each  forcing  mechanism 
on  the  mixed  layer  during  monsoon  onset. 
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TABLE  1 
PREVIOUS  STUDIES  OF  ATMOSPHERIC  FORCING 


Author  I Year 

Duing  and  Leetmaa  (1980) 

Krishnamurti(1981) 
Ramanadham  et  al.(1981) 


Rao  etal.  (1981) 


Pearce  and  Mohanty  (1984) 


Krishnamurti(1985) 


Gautier(1986) 


Study 

Discusses  heat  budget  for  the 
mixed  layer  between  April  and 
August;  Attributes  cooling  to 
coastal  upwelling  and  advection. 

Explores  role  of  strong  winds 
during  onset  in  1979  in  cooling 
Arabian  Sea.    Investigates  wind 
speed  and  wind  stress,  curl  of 
wind  stress  and  SST  fields,  1979. 

Investigates  energy  budget 
during  summer  monsoon,  1977. 
Analvzes  solar  radiation  flux, 
sensible  heat  flux,  evaporative 
heat  flux,  net  sea-surface  heat 

§ain,  wind  speed  and 
ST  fields. 

Examines  the  variability,  of 
mixed  layer  characteristics  during 
the  monsoon  in  1977. 

Examines  long  term  cooling 
of  the  Arabian  Sea  due  to  the 
summer  monsoon;  presents  wind 
speed,  net  moisture  flux  and 
net  heating  fields  for  1979, 
1980,  198fand  1982. 

An  overview  of  MONEX  that 
examines  the  wind  speed, 
wind  stress,  curl  of  wind 
stress,  SST  and  water  vapor 
flux  fields  during  1979. 

Describes  a  method  to 
derive  solar  radiation  at  the 
surface  and  describes  the 
evolution  of  this  field  during 
the  1979  onset. 
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II.  ANALYSES  OF  ATMOSPHERIC  FORCING  FIELDS 

A.       DISCUSSION  OF  PLAN 

This  study  will  address  several  possible  mechanisms  that  could  cause  the 
observed  decrease  in  the  SST  following  onset  of  the  summer  monsoon.  As  previously 
mentioned,  Krishnamurti  (1981)  suggested  five  possible  mechanisms.  Two  of  the  five 
mechanisms  will  be  discussed  and  analyzed,  along  with  two  other  possible  causes. 
These  cooling  mechanisms  include: 

(i)  Solar  radiative  flux; 

(ii)  Wind  speed-induced  mixing  and  entrainment; 

(iii)  Evaporative  heat  flux;  and 

(iv)  Sensible  heat  flux. 
As  Krishnamurti  (1981)  points  out,  the  northward  flux  of  heat  by  ocean  currents, 
as  estimated  by  Duing  and  Leetmaa  (1980),  is  small.  Furthermore,  this  advective  heat 
flux  could  not  respond  quickly  enough  to  account  for  the  rapid  cooling  that  occurs 
after  onset.  The  same  argument  can  be  applied  to  coastal  upwelling  of  cold  water. 
Therefore,  advective  heat  fluxes  are  not  considered  in  this  study.  Another  mechanism 
which  Krishnamurti  considers  important  for  cooling  is  wind-induced  upwelling.  Rao 
(1984)  analyzed  the  depth-time  variation  of  the  daily  averaged  isotherms  in  the  upper 
200  meters  of  the  Arabian  Sea  during  the  monsoon  onset  in  1979.  After  a  very  strong 
onset  event,  the  upward  velocity  of  the  isotherms  is  1  m/day.  The  heat  flux  at  the  base 
of  the  mixed  layer  is: 


Q(-h)  =  -cpp  ATdh/dt,  (2.1) 

where  Q(-h)  is  the  heat  exchange  at  depth  h  (the  mixed  layer  depth),  c  is  the  specific 
heat  of  seawater,  p  is  the  density  of  seawater,  AT  is  the  temperature  jump  across  the 
base  of  the  mixed  layer  and  dh/dt  is  the  vertical  velocity  of  the  interface,  which  is  given 
by: 


dh/dt  =  Wem  -  Wupw  ,  (2.2) 
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where  Went  is  the  entrainment  velocity  and  -W  w  is  the  upwelling  velocity.  For  a 
1°C  temperature  jump  and  an  upward  velocity  of  1  m/day,  the  heat  exchange  due  to 
upwelling  at  the  base  of  the  mixed  layer  is  approximately  50  W/m2.  This  is  a  small 
heat  flux  compared  to  other  fluxes  and  would  not  account  for  the  rapid  cooling  of  the 
mixed  layer.  The  1°C  temperature  jump  is  used  in  the  mixed  layer  prediction  part  of 
the  study  as  a  representative  value.  However,  if  the  temperature  jump  at  the  seasonal 
thermocline  was  greater,  vertical  advection  would  be  an  important  factor  in  the  cooling 
of  the  mixed  layer.  Although  vertical  advection  can  be  included  diagnostically  in  a 
one-dimensional  mixed  layer  prediction  model  (Muller  et  ai,  1984),  we  neglect  vertical 
and  horizontal  advection  in  this  study. 

B.       DATA  SOURCES 

This  study  utilizes  data  provided  by  the  NASA  Goddard  Space  Flight  Center  for 
May  and  June,  1979.  The  data  set  was  produced  from  the  NASA  four-dimensional 
data  assimilation  model  which  uses  the  previous  six-hour  prediction  field  as  a  first 
guess  for  the  current  analysis  field.  The  archived  fields  that  are  relevant  for  ocean 
prediction  studies  include:  sensible  heat  flux,  evaporative  heat  flux,  solar  radiative  flux, 
and  downward  fluxes  of  u-momentum  and  v-momentum.  The  assimilation  technique 
produces  output  four  times  per  day  at  a  grid  spacing  of  4°  latitude  x  5°  longitude. 

This  study  also  uses  similar  data  from  the  NOGAPS  model  of  FNOC  during 
May  and  June  of  1984  and  1985.  This  NOGAPS  output  is  also  produced  four  times 
daily,  with  analyses  at  00  and  12  GMT  and  six-hour  predictions  for  06  and  18  GMT. 
The  output  fields  include:  total  heat  flux,  sensible  heat  flux,  evaporative  heat  flux, 
solar  radiative  flux,  and  zonal  and  meridional  wind  speeds.  The  fields  are  on  a  73  x 
144  grid  with  2.5°  latitude  x  2.5°  longitude  spacing. 

These  data  are  analyzed  in  three  ways.  A  10-day  average  is  computed  for  the 
pre-onset  period  on  a  grid  centered  on  the  Arabian  Sea,  bounded  by  the  equator, 
25°N,  50°E  and  75°E  (Fig.  2.1).  The  pre-onset  period  for  1979  is  5  -  14  June.  The 
actual  onset  occurred  during  16-18  June.  The  pre-onset  period  for  1984  and  1985  is 
15  -  24  May  and  17  -  26  May,  respectively.  The  onset  occurred  during  26  -  28  May  for 
both  1984  and  1985.  A  10-day  average  is  also  computed  for  the  post-onset  periods  in 
all  three  years:  20  -  29  June  1979  and  3  -  12  June  in  both  1984  and  1985.  Standard 
deviations  from  these  average  values  for  the  pre-onset  and  post-onset  periods  are  also 
calculated  to  provide  a  measure  of  the  time  variability  of  the  atmospheric  forcing. 
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The  third  analysis  method  is  a  time  series  of  all  the  atmospheric  forcing  fields 
centered  on  the  onset  date.  The  time  series  are  examined  for  four  points  in  the 
Arabian  Sea  (Fig.  2.1)  that  are  located  well  away  from  any  currents.  The  northwest 
point  (NW)  is  at  15°N,  57.5°E;  the  northeast  point  (NE)  is  at  15°N,  67.5°E;  the 
southeast  point  (SE)  is  at  5°N,  67.5°E;  and  the  southwest  point  (SW)  is  at  5°N, 
57.5°E.  These  data  fields  will  be  compared  with  similar  data  from  other  studies  for  the 
purpose  of  characterizing  the  monsoon  onset  events  of  1979,  1984  and  1985  and  to 
verify  the  suitability  of  the  atmospheric  forcing  for  ocean  prediction  studies. 

C.       SURFACE  FLUXES 

The  heat  budget  equation  at  the  sea  surface  is: 

Qt  -  Qb  -  Qi  +  Qe  +  Qs .  <2-3) 

where  Qj  is  the  total  heat  flux  from  ocean  to  atmosphere;  Q^  is  the  net  flux  of  long- 
wave radiation;  Qj  is  the  downward  flux  of  solar  radiation;  Qe  is  the  upward  flux  of 
evaporative  (latent)  heat;  and  Qs  is  the  sensible  heat  flux.  The  Q^  term  is  not  explic- 
itly analyzed  in  this  study.  It  can  be  computed  as  a  residual  from  the  remaining  terms. 
The  long-wave  radiation  flux  actually  decreases  slightly  after  onset  because  of  increased 
cloud  cover,  but  the  magnitude  of  the  long-wave  radiation  flux  is  small  and  it  is  not  a 
major  factor  in  the  heat  budget.  These  heat  fluxes  are  summarized  in  the  following 
sections,  based  on  the  studies  in  Table  1. 
1.    Solar  Radiative  Flux 

This  flux  is  a  major  source  of  heating  prior  to  onset.  The  zenith  angle  of  the 
sun  is  very  low,  approximately  10°  on  May  15  at  10°N,  and  the  sky  is  generally  clear 
with  high  optical  transmissivity.  Gautier  (1986)  describes  how  the  solar  radiative  flux 
decreases  sharply  as  the  monsoon  onset  begins.  She  illustrates  the  formation  and 
northward  migration  of  the  Maximum  Cloudiness  Zone  (MCZ)  in  June,  1979.  The 
MCZ  clearly  follows  the  path  of  the  onset  vortex  that  is  traced  in  Krishnamurti  et  al. 
(1981)  and  Rao  (1984).  For  the  Arabian  Sea,  a  typical  value  of  daily  net  short  wave 
radiative  flux  before  onset  is  300  W/m2.  After  onset,  typical  flux  values  decrease  to 
about  250  W/m2  with  minimum  values  of  75  W/m2  under  the  MCZ.  Krishnamurti 
(1981)  indicates  that  the  rapid  cooling  of  the  Arabian  Sea  is  generally  not  attributed  to 
this  effect,  because  the  heat  capacity  of  the  ocean  causes  a  lag  time  of  one  to  two 
months  in  response  to  solar  radiation.   Nonetheless,  solar  radiative  flux  is  an  important 
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component  of  the  energy  balance  equation  as  will  be  shown  in  the  simulations  with  the 
Garwood  mixed-layer  model. 

The  solar  radiative  flux  values  are  available  each  six-hour  from  the  NASA  and 
NOGAPS  analyses.  This  time  interval  is  insufficient  to  indicate  the  true  maximum 
solar  radiation  as  a  function  of  time.  Therefore,  an  interpolation  scheme  (Gallacher, 
1979)  is  used  to  calculate  hourly  values  from  these  six-hour  values.  The  pre-onset  and 
post-onset  solar  radiative  fluxes  at  12  GMT  at  the  four  points  are  listed  in  Table  2  for 
1979,  1984  and  1985. 


TABLE  2 

MEAN  AND  STANDARD  DEVIATION  OF  SOLAR  RADIATIVE  FLUX, 

(W/M2) 

1979 

1984 

1985 

873 
(208) 

918 
(23) 

313 
(5) 

246 
(35) 

355 
(3) 

306 
(23) 

Pre- 
onset 

415 
(200) 

346 
(205) 

287 
(20) 

206 
(18) 

280 
(10) 

202 
(24) 

653 
(247) 

483 
(256) 

338 
(13) 

230 
(12) 

312 
(58) 

233 
(16) 

Post- 
onset 

693 
(288) 

527 
(282) 

340 
(10) 

197 
(3) 

311 
(24) 

210 
(10) 

Data    are    organized  by   year,    time   period,    location   and 
whether   mean   or    (standard   deviation).       For   exanrole,    in 
the    1979   pre-onset   period,    the   mean   fluxes    at   the   NW, 
NE,    SE    and    SW  points    are    873,    918,    346    and   415    W/m2, 

respecti1 

/e. 

Ly- 

The  solar  flux  values  are  generally  greater  to  the  north  during  the  northern 
summer  at  these  latitudes,  and  they  are  greater  to  the  west  where  12  GMT  lies  closer 
to  local  noon.  The  1979  (NASA  model)  solar  flux  values  appear  to  be  phase-shifted  as 
the  12  GMT  values  are  larger  than  the  06  GMT  values  (not  shown).  By  contrast,  the 
western  values  in  the  other  years  (FNOC  model)  are  larger  than  the  eastern  values. 
The  nature  and  reason  for  the  phase  difference  has  not  been  determined  in  this  study. 
Note  the  comparison  between  values  in  the  pre-onset  and  post-onset  periods.  The 
10-day  mean  values  at  the  northern  points  are  representative  of  cloudier  conditions 
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during  the  post-onset  period  than  are  the  southern  points,  where  little  or  no  decrease  in 
solar  flux  is  indicated.  The  1979  (NASA  model)  standard  deviations  are  significantly 
higher  than  are  the  1984  or  1985  values  (NOGAPS  model).  This  difference  is  attrib- 
uted primarily  to  the  difference  between  the  methods  used  to  calculate  solar  flux  in  the 
NASA  and  NOGAPS  models.  No  large  difference  is  evident  between  the  1984  and 
1985  data. 

The  means  and  standard  deviations  of  the  accumulated  daily  solar  radiative 
flux  (hourly  solar  radiation  flux  averaged  over  24  hours)  are  listed  in  Table  3.  These 
are  6-day  means  because  the  interpolation  of  the  solar  radiation  fluxes  did  not  include 
earlier  data.  The  accumulated  daily  solar  radiative  fluxes  are  computed  by  summing 
and  averaging  the  hourly  solar  radiation  flux  values  over  24  hours.  The  integrated 
solar  flux  values  for  all  three  years  are  very  comparable.  The  standard  deviations  are 
generally  10-15%  of  the  mean  values,  which  indicates  that  the  accumulated  solar  radi- 
ative fluxes  have  relatively  small  variability. 


TABLE  3 

AS  IN  TABLE  2, 

EXCEPT  FOR  AVG.  DAILY  SOLAR  RAD.  FLUX. 

1979 

1984 

1985 

Pre- 

111 

(12) 

124 

(16) 

165 
(19) 

193 
(20) 

158 
(18) 

233 
(22) 

onset 
Post- 

86 
(22) 

89 
(24) 

142 
(11) 

158 
(17) 

156 
(36) 

175 
(20) 

104 
(16) 

103 
(17) 

166 
(28) 

213 
(27) 

167 
(33) 

202 
(18) 

onset 

96 
(12) 

90 
(17) 

146 
(26) 

166 
(9) 

158 
(15) 

152 
(9) 

The  pre-onset  and  post-onset  12  GMT  instantaneous  solar  radiative  flux  fields 
averaged  over  10  days  from  17  to  26  May  and  3  to  12  June  1985  are  shown  in  Fig. 
2.2a  and  Fig.  2.2b,  respectively.    In  the  average,  there  appears  to  be  little  difference 
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between  the  solar  fluxes  in  the  two  periods.  As  shown  in  Fig.'  2.3,  the  time  series  of 
average  daily  solar  radiative  flux  at  the  SW  point  from  18  May  to  4  June  1985  indi- 
cates only  a  slight  decrease  of  daily  solar  radiative  flux.  This  is  also  true  at  the  other 
points  in  1984  and  1985.  In  1979,  when  an  onset  vortex  is  well  documented  in  various 
studies,  the  average  daily  solar  radiative  flux  at  the  northern  points  decreases  more 
than  in  the  other  years  (Table  3)  during  the  period  of  onset.  This  supports  the  findings 
of  Gamier  (1986).  It  is  possible  that  the  MCZ  in  1984  and  1985  simply  did  not  cross 
over  the  geographical  points  that  were  selected.  A  more  likely  possibility  is  that 
NOGAPS  does  not  properly  calculate  the  effect  of  cloudiness  on  solar  radiative  flux 
(see  Appendix  A).  In  any  case,  it  is  unlikely  that  the  increased  cloudiness  would 
decrease  the  solar  flux  over  'enough  area  and  a  sufficiently  long  time  to  cause  an 
appreciable  cooling  at  the  surface.  Therefore,  the  solar  radiative  flux  is  not  likely  to  be 
a  major  cause  of  the  observed  SST  decrease. 
2.    Wind  Mixing 

Several  studies  have  indicated  that  increased  wind  stress  should  induce  mixing 
and  entrainment  at  the  base  of  the  mixed  layer.  Ramanadham  et  al.  (1981)  illustrates 
the  wind  increase  at  four  points  in  the  Arabian  Sea  during  June  1977.  The  largest  wind 
increase  is  from  3  to  8  m/s  in  1  -  2  days  at  10.5°N,  66°E.  Krishnamurti  (1981)  shows 
an  even  stronger  increase  in  wind  speed  at  7°N,  66°E  during  the  onset  of  the  monsoon 
in  1979,  from  approximately  4  m/s  to  13  m/s  in  1  -  2  days. 

The  mixed  layer  warms  and  shallows  prior  to. the  monsoon  onset  as  a  result  of 
light  winds  and  strong  insolation.  Rao  et  al.  (1981)  examine  the  variability  of  the 
mixed  layer  in  the  area  of  the  central  Arabian  Sea  during  the  southwest  monsoon  in 
1977.  The  mixed  layer  depth  (MLD)  is  at  approximately  40  m  prior  to  onset  in  May 
and  June.  During  the  onset  phase,  the  MLD  increases  to  70  m  in  two  weeks  and 
deepens  further  as  the  monsoon  progresses.  The  mean  temperature  of  the  mixed  layer 
cools  by  1°C  within  one  week  after  onset.  Rao  et  al.  (1981)  attribute  this  deepening 
and  cooling  primarily  to  wind  mixing  and  also  to  the  convective  overturning  caused  by 
upward  heat  flux  at  the  surface.  It  is  the  hypothesis  of  this  study  that  the  increased 
wind  stress  at  the  sea-surface  mixes  the  surface  layer  and  entrains  cooler  water  from 
below  the  base  of  the  mixed  layer,  which  simultaneously  deepens  and  cools  the  layer. 

The  zonal  and  meridional  components  of  the  wind  are  extracted  from  the 
NOGAPS  analyses.  The  southwest  monsoonal  winds  generally  show  very  little  varia- 
tion  in  direction   after   onset,   so   the   variability   is   primarily  in  wind   speed.     The 
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calculation  of  wind  stress,  t,  is  described  in  Chapter  III.  Finally,  the  downward  fluxes 
of  u-momentum  (u'w')  and  v-momentum  (v'w')  are  computed  by  dividing  the  wind 
stress  components  by  the  air  density,  p.  These  fluxes  of  momentum  are  used  to 
compare  with  the  momentum  fluxes  which  have  been  extracted  directly  from  the 
NASA  analysis  for  1979.  Wind  stress  values  are  used  as  atmospheric  forcing  inputs  to 
the  Garwood  mixed  layer  prediction  model. 

Table  4  lists  means  and  standard  deviations  for  the  u-momentum  flux. 
Negative  values  indicate  a  downward  flux  of  westerly  momentum  since  w  is  negative  in 
the  downward  direction.  Notice  the  large  increase  in  westerly  momentum  flux  in  the 
post-onset  period  for  most  of  the  points.  Some  exceptions  are  the  momentum  fluxes  at 
the  SW  and  SE  points  in  1985  which  are  less  in  the  post-onset  period  than  in  the  pre- 
onset  period.  This  can  be  explained  by  noting  that  the  wind  speed  at  those  points  is 
stronger  after  onset  (see  Fig.  2.4a)  than  before  onset  (see  Fig.  2.4b),  but  the  westerly 
component  of  the  wind  is  less.  Notice  also  that  the  standard  deviations  are  generally 
larger  for  1979  than  for  1984  or  1985.  Since  the  data  for  1979  were  obtained  as  part  of 
FGGE  SOP  II,  it  is  likely  that  the  greater  number  of  observations  contributes  to  the 
larger  standard  deviations  in  1979,  whereas  the  NOGAPS  analyses  are  smoother  with 
smaller  standard  deviations.  However,  some  of  these  interannual  differences  may  be 
due  to  the  natural  variability  in  the  monsoon,  because  the  1979  monsoon  does  seem  to 
have  stronger  winds  overall. 

Table  5  lists  similar  calculations-  for  the  v-momentum  flux.  Negative  values 
indicate  a  downward  flux  of  southerly  momentum.  These  mean  values  are  generally 
less  than  the  corresponding  values  in  Table  4  because  the  monsoonal  winds  are 
predominantly  westerly,  especially  after  onset.  The  post-onset  periods  again  have 
significantly  higher  values  of  momentum  flux  than  do  the  pre-onset  periods.  In  partic- 
ular, the  v-momentum  fluxes  increase  sharply  in  1979.  The  onset  vortex  that  has 
already  been  discussed  is  responsible  for  these  effects.  It  is  obvious  from  these  two 
tables  that  the  monsoon  onset  is  accompanied  by  a  major  increase  in  surface  wind 
stress,  which  will  enhance  turbulent  mixing  in  the  ocean. 
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Pre- 
onset 


Post- 
onset 


TABLE  4 

AS  IN  TABLE  2.  EXCEPT  FOR  ATMOSPHERIC 
L -MOMENTUM  FLUX.  (M2  S*) 

19^9  1984  1985 


zSL-Q2        z0.02 
(0.09)       (0.04) 


z0.06 
(0.04) 


(0.20) 

zSL-11 

(0.06) 


0.00      ±0.02 

(0.02)       (0.05) 


^0.08 
(0.09) 


(0.02) 


^0.03 
(0.03) 


^0.0i 
(0.02) 

^0.05 
(0.03) 


(0.23) 

z0.16 
(0.14) 


zfi 

(0 


15 
10) 


zfi.ll 

(0.06) 


(0 


ii 
08) 


^0.07 
(0.03) 


zfl 

(0 


04 

05) 


z0-0i 
(0.02) 

H3-05 

(0.01) 


(0 


08 

05) 


nO.03      zfi -01 

(0.04)       (0.02) 


Two  representative  time  series  of  wind  speed  for  1985  are  illustrated  in  Fig.  2.5 
and  Fig.  2.6  for  the  SE  and  SW  points,  respectively.  At  the  SE  point,  the  wind  speed 
increases  by  2  m,  s  on  28  May  after  gradually  increasing  from  4  m  s  to  6.5  m  s  during 
the  previous  week.  At  the  SW  point,  the  wind  speed  steadily  increases  from  4  m  s  to 
13  ms  over  a  11-day  period.  The  axis  of  maximum  winds  is  aligned  closer  to  the  SW 
point  than  to  the  SE  point.  This  is  typical  since  the  low-level  jet  normally  extends 
southwest  to  northeast  along  the  Somali  coast  and  then  curves  eastward  across  the 
central  portion  of  the  Arabian  Sea.  Fig.  2.7  shows  the  evolution  of  u-momentum  flux 
at  the  NW  point  in  1979.  The  downward  flux  of  u-momentum  increases  abruptly  on 
17  June  and  a  similar  increase  in  v-momentum  flux  is  evident  in  Fig.  2.8.  Another 
interesting  feature  is  the  very  abrupt  change  in  sign  of  the  v-momentum  between  19 
and  21  June.  Based  on  the  path  of  the  onset  vortex  in  1979  (Krishnamurti  et  al.,  1981, 
and  Rao,  1984),  this  sudden  shift  in  winds  accompanies  passage  of  the  onset  vortex 
near  the  NW  point.  A  less  abrupt  shift  of  momentum  flux  is  seen  at  the  NE  and  SE 
points  a  few  days  earlier  as  the  onset  vortex  migrates  northward. 
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Pre- 
onset 


Post- 
onset 


TABLE  5 

AS  IN  TABLE  2,  EXCEPT  FOR  ATMOSPHERIC 
V-MOMENTUM  FLUX.  (M^S2) 


1979 


1984 


1985 


^0-05 
(0.08) 

(0.08) 


zQ.50 
(0.22) 


0 
(0 


00 
07) 


zO.fll 
(0.03) 

(0.03) 


0 
(0 

zfi 
(0 


oo 

01) 

02 
02) 


±0 
(0 

zSL 

(0 


04) 

02 
04) 


0 
(0 


00 
01) 


zo.oi 

(0.02) 


Z0.32 
(0.24) 


Z° 
(0 


02 
03) 


±0 
(0 


02 
04) 


^0 
(0 


05 
07) 


zP. 

(0 


02 
07) 


^0.09 
(0.04) 


^0.02 
(0.02) 


zl-08 
(0.04) 


±1 
(0 


01 
01) 


(0.02) 


^0.01 
(0.02) 

0.00 
(0.01) 


3.    Evaporative  Heat  Flux 

Ramanadham  et  al.  (1981)  show  that  the  evaporative  heat  flux  is  highly 
correlated  with  wind  speed,  with  increases  to  about  260  W/m2  as  the  monsoon  winds 
increase  and  values  of  about  150  W/m2  when  the  winds  subside.  Krishnamurti  (1981) 
lists  this  flux  as  possibly  being  an  important  factor  in  cooling  the  sea  surface. 

Table  6  is  a  compilation  of  the  10-day  mean  values  and  standard  deviations  of 
evaporative  heat  flux  for  the  four  points  during  the  pre-  and  post-onset  periods  of 
1979,  1984  and  1985.  The  mean  values  are  comparable  for  all  three  years,  but  the 
standard  deviations  of  the  1979  (NASA)  data  are  generally  significantly  higher  than  in 

1984  and  1985  (FNOC).  Also,  the  mean  values  for  the  northern  points  in  the  post- 
onset  period  are  approximately  twice  as  large  as  the  pre-onset  values.  The  fluxes  at  the 
southern  points  do  not  increase  significantly  (and  sometimes  decrease)  following  onset. 
This  follows  a  similar  trend  in  wind  speed  and  wind  stress  at  these  locations. 

Figs.  2.9a  and  2.9b  illustrate  the  distribution  of  the  evaporative  heat  flux  in 

1985  before  and  after  onset,  respectively.  Notice  the  spread  of  the  area  of  maximum 
flux  across  the  central  Arabian  Sea  after  onset.  Examples  of  the  evaporative  heat  flux 
increasing  with  time  at  the  SW  and  NW  points  during  1985  are  shown  in  Figs.  2.10 
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TABLE  6 

AS  IN  TABLE  2,  EXCEPT  FOR  EVAP. 

HEAT  FLUX.  (W/M2) 

1979 

1984 

1985 

Pre- 

132 
(80) 

91 
(71) 

112 
(1) 

222 
(18) 

161 
(18) 

132 

(1) 

onset 

151 
(49) 

120 
(58) 

381 
(17) 

339 
(9) 

342 
(50) 

277 
(3) 

Post- 

230 
(82) 

170 
(62)' 

295 
(18) 

350 
(28) 

379 
(9) 

321 
(33) 

onset 

189 

160 

382 

224 

405 

229 

(41) 

(84) 

(7) 

(0) 

(29) 

(25) 

and  2.11.    The  evaporative  heat  flux  increases  sharply  from  about  300  W/Mz  to  700 
W/M2  from  26  to  28  May  at  the  SW  point,  but  it  increases  by  only  about  150  W/Mz 
from  26  May  to  3  June  at  the  NW  point.    This  points  out  that  the  low-level  jet  is  a 
very  localized  feature  that  may  not  affect  all  locations  equally. 
4.    Sensible  Heal  Flux 

Ramanadham  et  al.  (1981)  analyze  the  sensible  heat  fluxes  over  the  eastern 
Arabian  Sea  during  June  1977.  The  magnitude  of  the  sensible  heat  flux  typically  is  on 
the  order  of  10  -  50  W/m2,  and  is  higher  at  the  time  of  onset  when  the  wind  speed  is 
greatest.  Krishnamurti  and  Ramanathan  (1982)  point  out  that  the  Bowen  ratio, 
expressed  as  sensible  heat  flux/evaporative  heat  flux,  is  typically  on  the  order  of  0.01  in 
the  tropics.  The  Bowen  ratio  in  the  study  by  Ramanadham  et  al.  (1981)  is  approxi- 
mately 0.1. 

The  10-day  average  Bowen  ratio  at  the  SW  point  at  12  GMT  is  120/600  =  0.2, 
which  indicates  the  relative  importance  of  the  evaporative  heat  flux  over  the  sensible 
heat  flux.  This  confirms  that  sensible  heat  flux  has  a  relatively  small  influence  on 
tropical  ocean  cooling  during  the  monsoon  onset. 

Table  7  lists  the  mean  values  and  standard  deviations  of  the  instantaneous  12 
GMT  sensible  heat  flux  values  for  the  NASA  and  NOGAPS  data.    Positive  values 
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indicate  sensible  heat  flux  from  the  ocean  to  the  atmosphere.  The  mean  values  of  flux 
are  all  positive  and  have  small  magnitudes.  They  do  not  show  a  large  increase  from 
the  pre-onset  period  to  the  post-onset  period.  The  NASA  data  for  1979  again  have 
higher  standard  deviations  than  do  the  NOGAPS  data  for  1984  and  1985. 


TABLE  7 

AS  IN  TABLE  2,  BUT  FOR  SENSIBLE  HEAT  FLUX.  (W/M2) 

1979 

1984 

1985 

17 
(16) 

7 
(8) 

20 
(10) 

57 
(9) 

43 
(7) 

20 
(4) 

Pre- 

onset 

28 

22 

58 

55 

57 

32 

(14) 

(10) 

(1) 

(0) 

(6) 

(3) 

27 
(22) 

13 
(22) 

38 
(21) 

55 
(9) 

17 
(7) 

44 
(11) 

Post- 

onset 

25 
(9) 

15 
(11) 

71 
(8) 

43 
(1) 

61 
(3) 

48 
(8) 

A  typical  time  series  of  sensible  heat  flux  is  given  in  Fig.  2.12  at  the  SW  point 
in  1985.  The  sensible  heat  flux  increases  significantly  during  the  onset  period.  The 
average  magnitude  of  the  flux  is  about  50  W/m2  before  onset  (27  May)  and  about  120 
W/m2  shortly  after  onset  (30  May).  This  increase  is  larger  than  Ramanadham  et  al. 
(1981)  found,  but  still  is  small  compared  to  the  evaporative  heat  flux  discussed  earlier. 
Sensible  heat  flux  is  therefore  not  expected  to  be  primarily  responsible  for  the  observed 
cooling. 

5.    Total  (Net)  Surface  Heat  Flux 

The  total  heat  flux  is  the  sum  of  the  three  heat  fluxes  discussed  above  plus  the 
downward  solar  radiative  flux  plus  the  net  flux  of  long-wave  radiation.  Since  the  total 
heat  flux  is  extracted  from  the  NOGAPS  atmospheric  prediction  model,  the  outgoing 
radiative  flux  can  be  calculated  from  the  remaining  terms.  This  study  does  not 
consider  the  long-wave  radiation  flux  for  analysis. 
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Ramanadham  et  al.  (1981)  indicate  the  total  heat  flux  from  ocean  to  atmos- 
phere increases  very  slightly  as  the  onset  begins,  and  decreases  strongly  as  the 
evaporative  and  sensible  heat  fluxes  decrease  after  onset.  However,  the  onset  in  1977 
appears  to  be  somewhat  anomalous  compared  to  the  values  in  this  study.  The 
NOGAPS  total  heat  flux  fields  for  1985  show  an  increase  in  the  daily  mean  heat  flux  of 
350  -  400  W/m2  from  the  pre-onset  to  the  post-onset  periods  at  the  SE  and  SW  points 
(see  Fig.  2.13  and  Fig.  2.14).  Positive  values  indicate  a  flux  of  heat  from  the  ocean  to 
the  atmosphere.  The  strong  diurnal  variation  of  the  total  heat  flux  is  obvious  in  these 
figures.  The  solar  radiation  during  the  day  forces  the  total  heat  flux  to  be  negative 
(into  the  ocean),  but  at  night  the  flux  is  to  the  atmosphere.  As  has  been  indicated,  the 
SW  point  is  affected  more  str&ngly  by  the  onset  and  the  total  heat  flux  increases  more 
at  this  point  than  at  the  SE  point.  These  flux  increases  are  consistent  with  the 
previous  analyses  of  the  component  heat  fluxes.  The  evaporative  heat  flux  is  the 
dominant  component  and  is  responsible  for  most  of  the  increase  in  total  heat  flux.  The 
two  northern  points  show  a  smaller  increase  of  approximately  150  -  200  W/m2  (see  Fig. 
2.15  and  Fig.  2.16).  The  corresponding  total  heat  flux  increases  in  1984  (not  shown) 
are  only  about  half  as  large,  which  indicates  a  more  gradual  or  less  intense  monsoon 
onset. 

D.      SEA-SURFACE  TEMPERATURE  FIELDS 

Several  studies  have  analyzed  SST  changes  during  the  onset  phase  of  the 
monsoon  in  the  Arabian  Sea.  The  SST  at  7°N,  66°E  decreases  from  30°C  to  27.5°C 
in  a  6-day  period  during  1979  as  the  wind  increases  from  4  to  14  m/s  (Krishnamurti, 
1981).  Ramanadham  et  al.  (1981)  indicate  that  the  SST  decreases  at  various  points  in 
the  Arabian  Sea  during  1977  are  on  the  order  of  0.5°C.  The  heat  fluxes  and  wind 
speeds  in  that  study  indicate  that  the  onset  was  not  dramatic,  which  explains  the 
smaller  decreases  in  SST  during  1977. 

The  NOGAPS  SST  data  are  produced  from  a  mix  of  observations,  climatology 
and  an  ocean  model  called  the  Thermal  Ocean  Prediction  System  (TOPS).  The  first- 
guess  field  for  the  SST  analysis  is  the  previous  12-hour  prediction  (Clancy  and  Pollack, 
1983).  In  the  absence  of  new  observations,  the  new  analysis  will  simply  be  the  first- 
guess  field  from  TOPS  with  a  small  correction  toward  climatology.  The  NASA  data 
set  available  to  us  does  not  include  the  SST  field  in  1979. 

The  time  series  of  SST  for  all  four  geographical  points  during  1984  and  1985  are 
given  in  Figs.  2.17  and  2.18,  respectively.    In  1984,  the  SST  at  the  NW  point  decreases 
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0.5°C  in  4  -  5  days  starting  1  June,  whereas  the  SST  at  the  NE  point  increases  and 
then  decreases  during  the  period.  A  stepwise  decrease  of  0.6°C  occurs  at  the  SE  point 
during  the  18-day  period.  The  SST  decreases  a  total  of  1.0°C  steadily  through  the 
entire  period  at  the  SW  point  in  1984.  In  1985,  the  NW  point  has  a  steady  decrease  of 
1.0°C  after  19  May.  A  constant  SST  is  found  at  the  NE  point  until  26  May  1985 
when  a  brief  increase  occurs  and  then  a  moderate  decrease  follows.  The  SW  point 
similarly  shows  a  modest  decrease  of  0.5°C  in  9  days  after  26  May  1985.  A  more 
dramatic  decrease  in  SST  of  nearly  1°C  in  8  -  9  days  starts  29  May  1985  at  the  SE 
point. 

The  SST  fields  used  by  the  NOGAPS  analyses  show  a  cooling  trend  at  the  NW, 
SW  and  SE  points  over  the  18-day  period  around  the  onset  in  1984  and  in  1985.  The 
time  and  rate  of  decrease  varies  from  point  to  point  and  year  to  year.  The  decrease  in 
the  sea-surface  temperature  is  not  as  great  in  1984  or  1985  as  Krishnamurti  (1981) 
reports  in  1979.  Only  the  decrease  of  1°C  in  8  -  9  days  at  the  SW  point  in  1985 
approaches  the  significant  cooling  event  of  1979.  No  onset  vortex  occurred  during 
1985.  The  monsoon  onset  of  1979  evidently  was  not  typical  in  strength  or  timing.  In 
both  1984  and  1985,  the  SST  decreases  during  onset  at  all  the  points  except  for  the  NE 
point.   Several  possibilities  exist  to  explain  this  anomaly: 


1.  The  NE  point  is  not  affected  by  the  onset  as  much.  However,  the  atmospheric 
forcing  analyses  do  not  support  this  statement.  The  onset  appears  to  be  more 
intense  than  at  the  NW  point  where  the  SST  decreases. 


The  cumulus  parameterization  and  solar  radiation  flux  calculations  are  inaccu 
rate.   This  may  be  ; 
results  and  does  no 


rate^  Thist  may  be  partially  correct,  but  again,  the  NW  point  should  show  similar 
t. 


3.  The  initial  mixed  layer  temperature  and  depth  are  more  resistant  to  onset 
cooling.  This  is  a  factor  that  cannot  be  verified  in  this  study  because  we  have  no 
verification  mixed  layer  data. 

4.  The  TOPS  model  does  not  parameterize  the  solar  absorption  in  the  mixed  layer 
correctly.  Again,  this  possibility  can  only  be  proven  if  verification  mixed  layer 
data  are  usedT 

The  reason  for  the  differences  in  the  NE  point  SST  field  remains  unanswered  in  this 

study,  but  the  solution  would  be  useful  for  mixed  layer  simulations  in  this  area. 

Which  of  the  atmospheric  forcing  fields  are  responsible  for  the  observed  SST 

cooling?    Much  of  the  following  discussion  will  focus  on  the  large  cooling  which  is 

observed  at  the  SE  point  in  1985.    From  the  atmospheric  forcing  analysis,  only  two 

mechanisms  appear  to  be  significant  factors  in  causing  the  SST  decrease:    evaporative 

heat  flux  and  wind  speed  or  stress.    Comparing  Figs.  2.10  and  2.6,  the  wind  speed  and 

evaporative  heat  flux  curves  increase  and  decrease  together  in  time  and  magnitude. 
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Increased  evaporative  heat  flux  causes  cooling  at  the  surface  and  increased  wind 
stress  causes  cooling  by  entrainment  at  the  base  of  the  layer.  The  Garwood  ocean 
mixed  layer  model  will  be  used  to  demonstrate  the  resulting  effect  on  SST  by  varying 
these  atmospheric  forcing  parameters. 
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Fig.  2.  1.   Locations  of  grid  points  in  the 

Arabian  Sea.   NOGAPS  anaysis  points:   A 
NASA  analysis  points:   ^  . 
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III.  LOCAL  FORCING  APPLIED  TO  AN  OCEAN  MIXED  LAYER  MODEL 

A.      THE  MODEL 

The  ocean  mixed  layer  model  (Garwood,  1977)  is  a  one-dimensional,  vertically 
integrated  model.  Consequently,  horizontal  and  vertical  advection  are  not  included  in 
the  model.  It  is  only  valid  for  predicting  mixed  layer  conditions  on  time  and  space 
scales  where  advection  can  be  neglected.  It  is  assumed  that  the  area  being  studied  is 
far  enough  from  the  Somali  current  and  the  time  scale  is  small  enough  (on  the  order  of 
5-10  days)  that  advection  is  not  a  significant  factor  for  cooling  the  mixed  layer.  The 
model  requires  four  sets  of  parameters  to  predict  the  response  of  the  mixed  layer:  the 
ocean  temperature  and  salinity  structure  in  the  upper  200  m;  net  surface  heat  flux  (Qn); 
magnitude  of  the  surface  wind  stress  (t);  and  the  solar  radiative  flux  (Qs).  The  three 
sets  of  atmospheric  forcing  data  are  obtained  directly  from  the  NOGAPS  model 
output.  The  net  heat  flux  is  the  sum  of  evaporative  heat  flux,  sensible  heat  flux,  solar 
radiative  flux  and  long-wave  radiative  flux.  Qn  is  positive  for  upward  heat  flux  from 
the  ocean  to  the  air.  As  was  discussed  in  Chapter  II,  the  evaporative  and  sensible  heat 
fluxes  both  increase  significantly  during  the  monsoon  onset.  Upward  fluxes  of  these 
quantities  contribute  directly  to  cooling  the  mixed  layer  at  the  surface  and  to  genera- 
tion of  convective  turbulence  that  cools  the  layer  by  entrainment  at  the  base  of  the 
layer. 

The  magnitude  of  the  surface  wind  stress  (t)  is  computed  from  the  NOGAPS 
wind  speeds  using  the  formulation  found  in  Gill  (1982): 
tx  =  CdpVu; 

Ty  =  CdPVv' 

T     =   (tx2  +   Ty2)l/2, 

where  tx  is  the  zonal  wind  stress;  Ty  is  the  meridional  wind  stress;  C^  is  the  drag  coef- 
ficient  for  the  ocean  surface;  p,  the  density  of  air,  is  set  to  1.22  kg/m  ,  V  is  the  wind 
speed;  and  u  and  v  are  the  zonal  and  meridional  components  of  the  wind  speed.  For 
speeds  at  or  below  6  m/s,  Cd  =  1.1x10  .  For  speeds  greater  than  6  m/s,  Cd  =  (0.61 
+  0.063V)xl0  .  Wind  stress  also  increases  significantly  during  the  monsoon  onset 
and  tends  to  deepen  the  mixed  layer  by  entrainment  of  water  from  below  the  mixed 
layer.  Ac  the  generation  of  turbulent  kinetic  energy  (TKE)  is  proportional  to  the  third 
power  of  the  wind  speed,  this  deepening  mechanism  is  highly  variable  even  after  onset. 
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Wind  stress  mixing  and  evaporative  heat  flux  are  obviously  closely  related.  However, 
the  sensitivity  study  can  indicate  whether  the  primary  heat  flux  takes  place  at  the 
surface  or  at  the  base  of  the  mixed  layer. 

As  described  in  Chapter  II,  the  solar  radiative  flux  (Qs)  does  not  decrease  as 
much  as  the  other  two  atmospheric  forcing  mechanisms  increase  during  onset.  In  the 
ocean  mixed  layer  model,  a  fixed  fraction  of  the  solar  radiation  is  absorbed  in  the  first 
meter  of  the  layer  and  the  remainder  of  the  absorption  occurs  exponentially  with  depth 
below  one  meter  (Stringer,  1983).  The  very  strong  insolation  before  onset  in  the 
Arabian  Sea  tends  to  produce  shallowing  and  warming  of  the  mixed  layer.  One 
purpose  of  this  study  is  to  determine  how  much  effect  the  increased  cloudiness  associ- 
ated with  onset  has  on  the  SST  and  mixed  layer  depth.  The  mixed  layer  model  uses 
the  wind  stress  and  surface  heat  flux  to  calculate  the  TKE  and  the  vertical  flux  of  TKE 
within  the  mixed  layer.  Garwood  (1977)  indicates  that  this  flux  is  the  major  source  of 
mixing  within  an  active  entrainment  zone  at  the  base  of  the  mixed  layer.  All  of  the 
data  are  input  on  an  hourly  basis.  Salinity  is  held  constant  with  depth  because  of  the 
lack  of  observational  data.  It  is  set  at  35.5  parts  per  thousand,  which  agrees  with 
surface  isohaline  distributions  found  in  Tchemia  (1981). 

Given  the  initial  structure  of  the  mixed  layer  and  the  atmospheric  forcing,  the 
model  predicts  the  changes  in  temperature  and  mixed  layer  depth.  Wind  stress 
produces  TKE  which  causes  mixing  and  entrainment.  However,  entrainment  offsets 
the  TKE  production  by  entraining  colder  water  with  less  TKE.  The  surface  buoyancy 
force  can  either  produce  TKE  (through  convective  overturning  initiated  by  surface 
cooling)  or  damp  TKE  (stabilization  due  to  surface  heating).  If  the  intensity  of  the 
TKE  and  surface  buoyancy  force  are  large  enough  to  overcome  the  buoyancy  force 
below  the  base  of  the  mixed  layer,  the  mixed  layer  will  deepen.  If  the  intensity  of  the 
TKE  or  the  surface  buoyancy  forces  diminish,  the  mixed  layer  will  retreat  to  a  shal- 
lower equilibrium  depth  on  the  order  of  the  Obukhov  depth.  This  depth  is  where 
surface  heating  of  the  mixed  layer  is  balanced  by  cooling  due  to  entrainment  or 
upwelling. 

The  atmospheric  forcing  terms  that  have  been  discussed  are  applied  to  this  ocean 
mixed  layer  model  to  investigate  three  questions: 

1.  How  well  does  the  Garwood  ocean  mixed  layer  model  predict  the  mixed  layer 
temperature  and  depth  using  the  atmospheric  forcing  terms  discussed  above  at 
low  latitudes  during  the  summer  months/ 

2.  Does  the  ocean  mixed  layer  model  predict  SST  changes  during  onset  that 
compare  favorably  with  the  SST  values  from  the  NOGAPs  analyses/ 
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3.   What  is  the  sensitivity  of  the  mixed  layer  predictions  to  each  of  the  three  atmos- 
pheric forcing  terms? 

The  monsoon  onsets  of  1984  and  1985  are  investigated  because  the  atmospheric 
forcing  terms  are  readily  obtainable  from  NOGAPS.  To  provide  the  initial  temperature 
fields  in  the  Indian  Ocean,  the  regional  Extended  Ocean  Thermal  Structure  (EOTS) 
analyses  for  1984  and  1985  were  interpolated  to  the  points  illustrated  in  Fig.  2.1.  The 
regional  EOTS  objective  analysis  system  (Clancy  and  Pollak,  1983)  integrates  all  avail- 
able Expendable  Bathythermograph  (XBT)  data  in  a  given  region  with  climatology.  It 
is  run  twice  a  week  with  all  XBT  observations  since  the  last  run.  EOTS  produces 
temperature  profiles  from  the  surface  to  400  m  at  17  fixed  levels  and  at  5  floating  levels 
which  are  centered  about  the  mixed  layer  depth.  The  horizontal  domain  is  a  63  x  63 
grid  with  approximately  50  km  grid  spacing  in  the  tropics  (Ken  Pollak,  FNOC, 
personal  communication).  When  XBT  observations  are  not  available,  the  analysis  is 
adjusted  toward  climatology.  A  significant  problem  is  the  paucity  of  observations  of 
both  atmospheric  and  oceanic  conditions  in  the  Arabian  Sea  region.  In  this  case,  very 
few  ocean  temperature  profiles  are  available  during  the  period  with  the  atmospheric 
forcing  data  from  NOGAPS. 

The  SST  from  the  EOTS  profiles  differ  in  magnitude  and  trend  relative  to  the 
NOGAPS  model  SST  analyses  described  above.  For  example,  the  NOGAPS  SST 
analysis  indicates  a  decrease  in  temperature  from  29.1  °C  to  28.3°C  at  the  SW  point  in 
1985.  The  corresponding  EOTS  SST  analysis  indicates  that  the  temperature  increases 
from  29.3  °C  to  29.5°C.  There  are  two  reasons  for  the  differences  in  SST.  First,  the 
two  data  sets  are  analyzed  on  two  different  grids.  The  NOGAPS  grid  spacing  is 
approximately  150  km  in  the  tropics  and  is  much  coarser  than  the  EOTS  grid. 
Secondly,  the  EOTS  data  indicate  no  apparent  response  to  the  onset  forcing  fields. 
Since  the  onset  occurs  at  a  different  time  each  year,  the  climatological  SST  analyses 
tend  to  average  out  the  response  to  onset  forcing.  The  EOTS  SST  analyses  appear  to 
be  based  mostly  on  climatology.  These  restrictions  prevent  a  verification  study  of  the 
ocean  thermal  structure  predictions.  Rather,  the  emphasis  will  be  on  the  relative 
sensitivity  of  the  mixed  layer  predictions  to  each  atmospheric  forcing  mechanism. 

Climatological  temperature  profiles  are  derived  from  Naval  Oceanographic  Office 
Reference  Publication  42C  (1982).  These  temperature  profiles  are  based  on  many  years 
of  XBT  observations  during  May  at  each  location.  The  average  SST  and  mixed  layer 
depth  are  used  to  construct  the  mixed  layer  profile.  Given  the  mixed  layer  temperature 
and  depth  and  the  temperature  at  200  m,  the  remainder  of  the  temperature  profile  is 
given  by: 
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T(z)  =  T0exp[Y(z-z0)],  (3.1) 

where  T(z)  is  the  temperature  at  depth  z;  Tq  is  the  mixed  layer  temperature;  Zq  is  the 
depth  of  the  mixed  layer,  and  y  =  ln[T(200)/T0]/[200-z0l.  The  temperature  at  1  m 
below  the  mixed  layer  depth  is  set  equal  to  Tq  -  1.  This  exponential  profile  approxi- 
mates the  climatological  profile  reasonably  well  (Fig.  3.1).  Thus,  the  initial 
temperature  profiles  are  determined  by  the  SST  (which  is  also  the  mixed  layer 
temperature),  the  mixed  layer  depth  and  the  temperature  at  200  m,  as  listed  in  Table  8. 


TABLE  8 

Climatological  SST, 
temperature  at  200m 

mixed  layer 
at  the  four 

depth  and 
points. 

SE 

SW      NW 

NE 

SST  (*C)          29.9 

29.3    29.0 

31.9 

M.L.  Depth  (m)     40 

60      30 

30 

T200  (°C)         14.2 

14.8     15.5 

17.2 

The  shallower  mixed  layer  at  the  northern  two  points  relative  to  the  two  southern 
points  is  probably  because  the  axis  of  the  developing  low-level  jet  is  closer  to  the 
southern  points  just  prior  to  onset.  The  SST  at  the  NE  point  is  suspiciously  high 
compared  to  the  climatological  isotherm  analyses  of  Tchernia  (1981)  and  the  Monthly 
Oceanography  Summary  for  May  1985.  In  any  case,  the  initial  SST  has  only  a  limited 
effect  on  the  mixed  layer  model.  The  simulated  mixed  layer  responds  quickly  to  the 
atmospheric  forcing  and  an  equilibrium  thermal  structure  develops  within  4-5  days. 
The  lower  SST  at  the  NW  point  is  probably  due  to  the  influence  of  the  Somali  Current 
as  it  spreads  over  the  northern  Arabian  Sea.  The  actual  SST  at  the  NE  point  is  prob- 
ably somewhat  higher  than  the  SST  at  the  southern  points  because  of  the  lighter  winds. 
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B.       MIXED  LAYER  PREDICTIONS  DERIVED  FROM  THE  ATMOSPHERIC 
FORCING 

The  atmospheric  forcing  fields  are  applied  to  the  ocean  mixed  layer  model  in  two 
phases.  First,  the  FNOC  forcing  fields  during  1984  and  1985  are  applied  at  all  four 
points  to  determine  the  mixed  layer  development  and  to  compare  the  predicted  SST  to 
the  NOGAPS  analyzed  SST.  A  similar  test  of  the  NASA  atmospheric  forcing  during 
1979  is  not  done  as  no  verification  data  on  the  evolution  of  the  SST  are  available  to  us. 
Second,  pre-onset  mean  atmospheric  forcing  values  are  substituted  for  the  actual  data 
to  determine  the  evolution  of  the  mixed  layer  temperature  and  depth  in  the  absence  of 
an  onset  event.  Each  forcing  mechanism  is  then  examined  by  selectively  substituting 
the  actual  forcing  values  while  keeping  the  other  two  forcing  mechanisms  at  their  pre- 
onset  mean  values.  The  second  phase  is  illustrated  for  the  SE  point  only  during  the 
period  22  May  to  6  June  1985.  This  point  is  chosen  because  it  is  an  area  with  wide 
variations  between  the  pre-onset  and  the  post-onset  atmospheric  forcing  fields.  The 
integrations  are  begun  on  16  May  to  ensure  that  equilibrium  has  been  achieved  with 
five  days  of  pre-onset  forcing  before  the  output  is  examined.  The  onset  of  the 
monsoon  was  about  26  -  27  May  in  1984  and  1985. 

The  evolution  of  the  mixed  layer  temperature  and  depth  during  the  1985  onset  is 
illustrated  in  Figs.  3.2  and  3.3,  respectively.  Daily  maximum  mixed  layer  depths  are 
presented,  along  with  the  corresponding  mixed  layer  temperature  at  that  time  (usually 
early  morning),  because  the  diurnal  changes  predicted  by  the  mixed  layer  model  cannot 
be  verified.  At  the  SE  point,  the  pre-onset  SST  is  29.4°C.  and  the  mixed  layer  depth 
is  45  m.  Beginning  on  26  May,  the  SST  decreases  steadily  to  about  28.3°C  and  the 
mixed  layer  deepens  to  57  m.  As  seen  in  Figs  2.10  and  2.16,  this  is  approximately  the 
date  that  the  net  heat  flux  and  wind  speed  increase  in  magnitude. 

At  the  SW  point,  the  pre-onset  SST  is  29.0°C  and  the  mixed  layer  depth  is  62  m. 
Following  monsoon  onset,  the  mixed  layer  cools  by  1.5°C  and  deepens  by  21  m  in  nine 
days.  This  large  cooling  and  deepening  at  the  SW  point  is  even  more  significant  since 
the  initial  mixed  layer  depth  was  largest  at  this  point. 

Prior  to  onset,  the  mixed  layer  at  the  NW  point  warms  and  shallows  very  rapidly 
in  response  to  atmospheric  forcing.  The  average  mixed  layer  temperature  is  32°C  and 
the  mixed  layer  depth  is  3-8  m.  Formation  of  such  a  warm,  shallow  layer  at  the 
surface  is  an  excessive  response  in  the  model  to  the  strong  insolation  and  light  winds 
that  exist  prior  to  onset.  If  the  SST  was  really  this  high,  the  surface  heat  fluxes  would 
be  much  higher  than  the  NOGAPS  atmospheric  forcing  fields,  which  are  based  on 
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actual  SST  values.  That  is,  a  fully  coupled  ocean-atmosphere  model  in  which  the 
surface  fluxes  depend  on  the  predicted  SST  would  be  more  realistic  than  the  procedure 
used  here  of  prescribing  externally  the  actual  atmospheric  forcing  fields. 

There  are  several  possible  explanations  for  the  excessive  warming  and  shallowing 
that  is  simulated  in  the  mixed  layer  model.  The  NOGAPS  analysis  uses  an  algorithm 
to  indirectly  estimate  the  cloud  cover  based  on  the  humidity  fields.  The  algorithm 
could  result  in  too  high  values  for  Qs  which  would  cause  too  much  warming  and  shal- 
lowing. However,  the  NOGAPS  solar  fluxes  described  above  appear  to  be  too  small 
relative  to  Gautier  (1986).  The  mixed  layer  model  may  not  use  the  proper  distribution 
of  solar  radiation  absorption  in  the  mixed  layer.  If  too  large  a  fraction  of  the  Q  is 
assumed  to  be  absorbed  in  the  first  meter,  it  would  also  cause  excessive  warming  and 
shallowing.  There  are  usually  very  few  wind  observations  in  this  area  for  use  in  the 
NOGAPS  analysis.  Low  wind  speeds  from  NOGAPS  could  also  result  in  a  shallow 
warm  layer.  As  there  may  be  a  combination  of  causes  in  this  case,  the  answer  is 
unclear.  Stringer  (1983)  reported  similar  predictions  of  warm  and  shallow  layers  during 
the  spring  and  summer  months  in  the  North  Pacific  Ocean. 

Beginning  on  25  May,  the  mixed  layer  at  the  NW  point  immediately  starts 
cooling  and  deepening.  Although  the  net  temperature  difference  between  22  May  and 
6  June  is  nearly  zero,  the  layer  cools  by  2.5°C  and  deepens  by  20  m  between  25  May 
and  6  June.  During  this  period  the  fictitious  warm  and  shallow  layer  near  the  surface 
is  eliminated. 

The  NE  point  experiences  even  more  excessive  warming  prior  to  onset  and 
similar  cooling  and  deepening  after  onset.  The  net  mixed  layer  temperature  and  depth 
change  (positive  values  represent  cooling  and  deepening)  as  well  as  the  maximum 
mixed  layer  temperature  and  depth  change  within  the  period  are  listed  in  Table  9  for  all 
four  points.   The  results  for  the  1984  onset  period  show  the  same  tendencies  as  in  1985. 

Simpson  and  Dickey  (1981)  examined  two  classes  of  parameterization  of  solar 
irradiance  in  various  mixed  layer  models.  The  first  class  used  a  single  attenuation 
length  and  the  other  class  of  models  used  two  attenuation  lengths.  The  models  with 
two  attenuation  lengths  were  more  sensitive  to  the  surface  fluxes  and  absorbed  more 
solar  radiation  in  the  upper  few  meters.  The  differences  between  the  two  classes  of 
models  were  greater  for  light  wind  conditions,  which  is  the  case  during  the  pre-onset 
period  at  the  northern  points.  Martin  (1985)  found  that  the  turbidity  also  affected  the 
absorptive  properties  of  the  mixed  layer.    These  studies  indicate  that  the  mixed  layer 
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model  must  be  adjusted  to  the  conditions  of  the  atmosphere  and  ocean  if  it  is  to  simu- 
late accurately  the  mixed  layer  changes.  The  Garwood  mixed  layer  model  uses  two 
attenuation  lengths  (one  of  which  is  very  short),  and  is  not  adjusted  to  account  for 
possible  variations  in  the  turbidity  of  the  water.  The  model  is  tuned  for  midlatkude 
solar  fluxes,  wind  speeds  and  turbidity,  and  it  is  not  able  to  simulate  accurately  the 
pre-onset  mixed  layer  changes. 


TABLE    9 

Initial    SST    and  mixed   layer   depth   and   resulting   cooling 
(    C)    and   deepening   (m)    predicted  by   the   mixed 
layer   model    using    actual    forcing  values    at    the    four 
points    in    1984    and    1985.       Net   value    is    the   difference 
between   the    final    and   initial   value.    Max   value    is    the 
difference   between   the   max.     and  min.    values    in   the   period. 

YEAR 

PT 

PRE-ONSET 
SST        MLD 

COOLING 
NET        MAX 

DEEPENING 
NET      MAX 

1984 

SE 

29.5 

42 

0.4 

0.4 

6 

37 

sw 

29.0 

63 

0.9 

0.9 

11 

11 

NW 

30.0 

3 

0.4 

0.3 

27 

28 

NE 

32.4 

3 

0.2 

1.5 

28 

28 

1985 

SE 

29.0 

45 

1.1 

1.1 

11 

11 

SW 

29.0 

62 

1.5 

1.5 

20 

27 

NW 

31.7 

3 

0.0 

2.5 

19 

19 

NE 

36.7 

2 

1.8 

2.8 

10 

10 

Table  10  lists  the  values  for  SST  for  both  the  ocean  mixed  layer  model  and  the 
NOGAPS  analysis  for  all  four  points  before  and  after  onset  in  1985.  Only  at  the  SE 
point  do  the  NOGAPS  and  mixed  layer  model  results  agree.  Despite  using  a  similar 
temperature  profile  to  initialize  the  mixed  layer  model,  the  atmospheric  forcing  mecha- 
nisms at  the  northern  points  cause  the  mixed  layer  to  warm  and  shallow  rapidly  during 
the  pre-onset  period  so  that  the  NOGAPS  SST  and  mixed  layer  temperatures  are  quite 
different  when  onset  begins.    For  these  reasons,  this  study  will  analyze  the  sensitivity  of 
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TABLE    10 

Pre-onset    and   post-onset    SST    ( 
layer   model    and  NOGAPS    at    the 

°C)    for   the   mixed 
four  points    in    1985. 

Pre-onset    SST 
ML    Model             NOGAPS 

Post-onset    SST 
ML    Model           NOGAPS 

SE                   29.4 

29.4 

28.3               28.4 

SW                    29.0 

29.2 

27.5               28.4 

NW                    31.7 

28.5 

31.7               27.7 

NE                    36.7 

28.5 

35.0               28.3 

the  mixed  layer  predictions  to  the  forcing  terms  rather  than  emphasizing  the  compar- 
ison between  the  mixed  layer  model  SST  and  the  NOGAPS   SST. 

In  Phase  2a,  constant  values  for  Qn,  T  and  Qs  representing  pre-onset  conditions 
are  provided  to  the  model.  The  purpose  of  this  procedure  is  to  establish  a  baseline  for 
comparing  the  relative  effects  of  each  of  the  three  atmospheric  forcing  terms.  The 
primary  requirement  is  that  the  mixed  layer  be  in  a  relatively  steady  state  while  using 
reasonable  values  for  the  forcing  terms.  If  the  actual  pre-onset  values  are  used,  the 
mixed  layer  would  gradually  warm  since  the  strong  insolation  would  not  be  balanced 
by  the  net  surface  heat  losses  or  by  mixing  due  to  wind  stress.  These  pre-onset  values 
are  based  on  the  results  of  the  atmospheric  forcing  analysis  presented  in  Chapter  II. 
However,  they  have  been  empirically  adjusted  to  maintain  a  nearly  constant  mixed 
layer  temperature  and  depth  throughout  the  pre-onset  period.  The  values  are:  Qn  = 
135  W/m2;  T  =  0.02  N/m2,  and  Qs  =  150  W/m2.  Therefore,  the  net  surface  heat  flux 
is  approximately  50%  larger  than  the  actual  pre-onset  value  and  the  other  forcing 
terms  are  representative  of  the  pre-onset  values.  The  resulting  mixed  layer  characteris- 
tics are  listed  in  Table  1 1  . 

With  this  baseline  pre-onset  structure  established,  the  sensitivity  of  the  mixed 
layer  predictions  to  each  forcing  term  is  tested  in  Phases  2b-2d  (  Figs.  3.4  and  3.5).  In 
Phase  2b,  the  effect  of  Qn  is  tested  by  inserting  actual  values,  while  T  and  Qs  are  held 


53 


constant  at  the  pre-onset  values.  As  seen  in  Table  11,  the  net  surface  heat  flux  has  a 
very  significant  impact  on  the  cooling  and  deepening  of  the  layer.  Comparing  values  in 
Table  11,  the  actual  SST  decrease  (Fig.  3.4)  is  1.1  °C  and  the  Qn  term  alone  accounts 
for  1.1  °C  cooling.  Qn  also  accounts  for  almost  half  of  the  mixed  layer  deepening  that 
occurs  when  complete  atmospheric  forcing  is  specified. 

In  Phase  2c,  the  actual  values  of  surface  wind  stress  are  inserted  while  Qn  and  Qs 
are  held  constant  at  the  pre-onset  values.  The  increase  in  surface  wind  stress  after 
onset  has  less  effect  than  did  the  Qn  forcing.  The  SST  decreases  by  0.2°C  and  the 
mixed  layer  deepens  by  8  m  when  actual  values  for  T  are  used  in  the  model.  Even 
though  the  mixed  layer  deepens  nearly  as  much  as  it  does  with  the  Qn  forcing,  the 
cooling  is  much  less  because-  the  larger  direct  surface  cooling  after  onset  is  not 
included.  Mechanical  mixing  and  subsequent  entrainment  alone  do  not  appear  to  be 
responsible  for  most  of  the  cooling  of  the  mixed  layer.  However,  at  least  half  of  the 
deepening  can  be  ascribed  to  wind  mixing. 

When  actual  values  for  Qs  are  used  (Phase  2d),  the  mixed  layer  warms  0.2°C  and 
deepens  3  m.  One  would  expect  Qs  to  decrease  during  onset  (Gautier,  1986),  but  the 
Qs  values  from  NOGAPS  do  not  decrease.  Thus,  the  question  of  the  relative  effect  on 
the  ocean  mixed  layer  of  the  insolation  decrease  during  monsoon  onset  can  not  be 
addressed  with  these  NOGAPS  values. 


TABLE 

11 

As    in  Table    9,    except    for    1985    at   the    SE   pt.     only: 
Phase   2a.    using   all   pre-onset    forcing  values;    phase   2b, 
using   actual    Q      forcing;    phase   2c ,    using   actual 
T   forcing;    phasfe   2d,    using   actual   Q_    forcing. 

PHASE 

YEAR 

PT 

PRE- 
SST 

DNSET 
MLD 

COOLING           DEEPENING 
NET      MAX        NET      MAX 

2a 

1985 

SE 

30.0 

41 

0.0 

0.0           1 

1 

2b 

1985 

SE 

29.6 

42 

1.1 

1.1           9 

9 

2c 

1985 

SE 

29.9 

42 

0.2 

0.2          8 

8 

2d 

1985 

SE 

30.0 

33 

-0.2 

-0.2           3 

10 
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In  summary,  the  results  of  Phase  1  indicate  the  ocean  mixed  layer  model  appears 
to  simulate  successfully  mixed  layer  changes  during  the  onset  at  the  SE  and  SW  points. 
The  predictions  using  the  atmospheric  forcing  at  the  NW  and  NE  points  produce 
unrealistic  SST  values  when  compared  to  other  climatological  data.  The  fault  may  be 
with  the  mixed  layer  model  or  with  the  NOGAPS  forcing  fields  or  both.  An  exact 
answer  can  only  be  found  if  the  forcing  fields  and  the  resulting  mixed  layer  profiles  are 
compared  to  true  verification  data  sets.  In  Phase  2,  the  net  surface  heat  flux  term  has 
the  largest  influence  in  cooling  the  mixed  layer  during  onset,  but  the  net  surface  heat 
flux  and  the  wind  stress  contribute  almost  equally  to  mixed  layer  deepening.  Although 
the  net  surface  heat  flux  and  wind  stress  are  not  independent  and  cannot  be  totally 
separated,  the  model  results  suggest  that  the  heat  flux  at  the  surface  is  greater  than  at 
the  base  of  the  mixed  layer.  That  is,  while  entrainment  causes  significant  deepening, 
the  water  that  is  entrained  does  not  cool  the  layer  significantly. 
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Fig.  3. 1.   Exponential  ( solid)  approximation  of 
the  climatological  (dashed)  temperature 
profile  at  the  SE  point  during  May. 
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Fig.     3. 2.       Mixed    layer   temperature    at   time   of  maximum 
mixed    laver   depth   each   day    at    all   points,     1985. 
Solid   -    SE   point;    Dash   -    SW   point; 
Chain-dot   -    NW  point;    Chain-dash   -   NE   point. 
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Fig.  3.  3.   Daily  maximum  mixed  layer  depth  with  1985  forcing 

values  applied.   Same  line  pattern  as  in  Fig.  3.2. 
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Fig.     3. 4.       Mixed    layer   temp,     at   time   of   max.    mixed    layer 
depth   each   day  with  various    1985    post-onset    forc- 
ing  applied.       Solid:     pre-onset    forcing    (Phase    2a); 
dash:     Q_    forcing   (Phase    2b);    chain-dot:     t    forc- 
ing  ( Phise    2c);    chain-dash:     Qs    forcing    (phase    2d). 
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3.5.       Daily  max.    mixed    layer   depths   with  various    1985 
post-onset    forcing    applied.     Same   key   as    in   Fig.     3.  4. 
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IV.  CONCLUSIONS 

A.  BACKGROUND 

The  purpose  of  this  study  has  been  to  determine  the  most  significant  causes  for 
the  rapid  sea-surface  temperature  decrease  in  the  Arabian  Sea  during  the  onset  of  the 
southwest  monsoon.  The  atmospheric  forcing  of  the  central  Arabian  Sea  at  the  time  of 
onset  has  been  examined.  Atmospheric  forcing  data  for  1979  are  available  from  the 
NASA  four-dimensional  data  assimilation  model,  and  data  for  1984  and  1985  are  taken 
from  the  NOGAPS.  The  forcing  fields  have  been  analyzed  at  four  points  in  the 
Arabian  Sea  to  determine  the  pre-onset  and  post-onset  values.  The  NOGAPS  forcing 
fields  are  chosen  to  apply  to  the  Garwood  ocean  mixed  layer  model  because  no  SST 
fields  for  1979  are  available  for  verification.  The  model  is  initialized  with  approxima- 
tions to  the  climatological  temperature  and  salinity  profiles.  The  model  is  then  run  for 
the  period  covering  the  monsoon  onset.  Phase  1  of  this  procedure  is  the  application  of 
the  actual  onset  forcing  parameters  to  determine  how  well  the  ocean  mixed  layer  could 
be  simulated  given  these  data.  Phase  2  involves  the  application  of  fixed  pre-onset 
forcing  values  to  the  model.  The  actual  onset  values  of  each  forcing  parameter  are 
then  substituted  one  at  a  time  to  test  the  sensitivity  of  the  mixed  layer  predictions  to 
each  forcing  term.  The  results  are  also  examined  in  terms  of  differences  among  the 
four  points. 

B.  ATMOSPHERIC  FORCING  ANALYSIS 

The  atmospheric  forcing  fields  are  evaporative  heat  flux,  wind  speed,  solar  radia- 
tive flux  and  sensible  heat  fiux.  The  long-wave  radiative  flux  is  not  examined  because 
it  does  not  change  appreciably  with  the  onset  and  is  not  expected  to  contribute  to  the 
unusual  cooling  of  the  mixed  layer.  It  is  important  to  realize  that  these  forcing  terms 
are  not  observed  quantities,  but  are  the  output  of  a  four-dimensional  data  assimilation 
system.  The  NOGAPS  analysis  has  relatively  poor  observational  input  in  this  region 
and  has  a  coarse  grid  resolution.  Therefore,  the  atmospheric  forcing  fields  must  be 
utilized  with  caution. 

The  solar  radiative  flux  appears  to  be  very  dependent  upon  the  algorithm  used  to 
estimate  it  in  the  atmospheric  model.  Although  the  solar  radiative  flux  from  the 
NASA  analysis  is  phase-shifted  in  time,  the  total  daily  heat  fluxes  are  comparable  to 
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the  NOGAPS  solar  radiative  flux.  Whereas  Gautier  (1986)  found  a  large  decrease  in 
insolation  over  much  of  the  Arabian  Sea  after  the  1979  onset,  neither  the  NASA  anal- 
ysis for  1979  nor  the  NOGAPS  analyses  for  1984  and  1985  indicate  a  significant 
decrease  in  solar  radiative  flux  at  any  of  the  points.  Therefore,  the  calculations 
involving  solar  radiative  flux  as  a  possible  cause  of  the  mixed  layer  cooling  are  incon- 
clusive. The  method  for  calculating  this  parameter  in  the  NOGAPS  analysis  should  be 
examined  as  the  potential  cause  of  erroneous  ocean  predictions  in  the  tropics. 

The  evaporative  heat  flux  increases  dramatically  during  onset.  It  is  the  most 
significant  upward  heat  flux  after  onset  and  it  dominates  the  total  (net)  surface  heat 
flux.  Much  of  the  variation  in  evaporative  heat  flux  is  associated  with  the  changes  in 
wind  speed. 

The  sensible  heat  flux  more  than  doubles  in  magnitude  after  the  onset,  but  it  is 
still  significantly  smaller  than  evaporative  heat  flux  and  solar  radiative  flux.  It  is  not 
believed  that  the  sensible  heat  flux  contributes  significantly  to  the  mixed  layer  cooling 
during  onset. 

The  final  atmospheric  forcing  mechanism  is  the  wind  stress.  It  is  appears  to  be 
the  most  immediate  indicator  of  the  monsoon  onset  over  the  Arabian  Sea.  The  wind 
speed  dramatically  increases  as  the  monsoon  begins  in  all  three  years.  This  increase  is 
most  obvious  near  the  core  of  the  low-level  Somali  jet,  or  to  the  south  of  the  anticy- 
clonic  onset  vortex  as  in  1979.  However,  it  does  increase  to  some  extent  over  the 
entire  Arabian  Sea.  The  increased  wind  speed  has  two  effects  on  the  ocean  mixed 
layer.  First,  the  increased  wind  stress  generates  more  turbulence  in  the  mixed  layer, 
which  will  entrain  cooler  water  from  below.  Second,  the  higher  wind  speed  increases 
the  evaporative  and  sensible  heat  fluxes.  In  particular,  the  wind  speed  and  evaporative 
heat  flux  varied  in  the  same  manner  during  the  onset  period.  For  this  reason,  it  is 
difficult  to  distinguish  the  relative  importance  of  evaporative  heat  flux  and  wind  speed. 

C.       EXAMINATION    OF    ATMOSPHERIC    FORCING    WITH    THE    MIXED 
LAYER  MODEL 

Hourly  values  of  the  net  surface  heat  flux,  surface  wind  stress  and  solar  radiative 

flux  are  derived  from  the  NOGAPS  model  output  for  1984  and  1985  and  used  as  input 

to  a  one-dimensional,  ocean  mixed  layer  model.    The  net  surface  heat  flux  is  the 

combination  of  sensible  heat  flux,  solar  radiative  flux,  long-wave  radiative  flux  and 

evaporative   heat   flux.     Climatological   temperature   profiles   and   isohaline    salinity 

profiles  are  used  for  the  initial  mixed  layer  characteristics  as  accurate  profiles  are  not 
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available  from  the  operational  analyses.  The  model  is  initialized  on  16  May  of  both 
years,  which  is  approximately  11  days  before  onset,  and  is  run  until  6  June,  which  is 
approximately  10  days  after  onset. 

At  the  southern  points,  the  mixed  layer  model  is  able  to  predict  SST  changes 
which  are  close  to  the  NOGAPS  SST  changes.  The  SST  decrease  starts  immediately 
after  the  first  indications  of  onset  based  on  the  atmospheric  forcing  fields.  However, 
the  pre-onset  mixed  layer  temperatures  at  the  northern  points  are  unrealistically  high 
after  five  days  of  integration.  The  solar  radiative  flux  from  NOGAPS  may  be  inaccu- 
rate or  the  mixed  layer  model  may  not  integrate  the  forcing  terms  properly  under 
conditions  of  strong  insolation  and  light  winds.  More  observations  are  needed  to 
resolve  this  problem,  which'  severely  limits  the  usefulness  of  ocean  mixed  layer 
predictions  in  this  northern  area. 

Based  on  a  sensitivity  study,  the  net  surface  heat  flux  is  the  most  significant 
mechanism  contributing  to  mixed  layer  cooling.  When  the  other  forcing  terms  are 
fixed  at  pre-onset  values,  the  post-onset  values  of  the  net  heat  flux  account  for  as 
much  cooling  (1.1°C)  as  when  all  the  forcing  terms  are  set  at  post-onset  values.  The 
net  heat  flux  and  the  wind  stress  have  nearly  equal  contributions  to  the  predicted 
deepening  after  onset.  Increased  wind  stress  alone  forces  some  deepening  but  only 
modest  cooling  of  the  mixed  layer  (0.2°C).  Since  evaporative  heat  flux  is  a  major  part 
of  the  net  heat  flux,  the  latter  is  necessarily  tied  to  changes  in  wind  speed  (stress). 
However,  the  results  do  imply  that  heat  loss  from  the  sea  surface  is  greater  than  the 
downward  heat  flux  at  the  base  of  the  mixed  layer  due  to  entrainment  mixing.  This 
suggests  that  the  original  hypothesis  of  this  study  is  must  be  modified.  Both  the  wind- 
induced  mixing  and  the  evaporative  heat  flux  must  be  considered  to  explain  the  mixed 
layer  cooling  and  deepening  in  the  Arabian  Sea  during  the  monsoon  onset. 

In  summary,  this  study  has  only  been  partially  successful  in  determining  the 
oceanic  effects  of  changes  in  the  atmospheric  forcing  terms  during  onset  of  the  south- 
west monsoon.  The  conclusion  regarding  the  importance  of  net  surface  heat  flux  must 
be  regarded  as  tentative  because  the  algorithm  used  by  NOGAPS  and  NASA  to  calcu- 
late the  solar  radiative  flux  may  be  unsuitable  under  these  southwest  monsoon 
conditions.  Given  accurate  atmospheric  forcing,  accurate  ocean  thermal  structure 
observations  for  verification  and  a  properly  adjusted  ocean  mixed  layer  model,  the 
study  suggests  that  the  ocean  mixed  layer  changes  during  monsoon  onset  can  be 
accurately  predicted. 
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This  analysis  is  a  useful  first  step  in  understanding  the  cause  of  the  cooling  of  the 
Arabian  Sea  during  the  monsoon  onset.  Future  studies  should  use  atmospheric  forcing 
based  more  on  observations  rather  than  on  model  predictions.  The  data  from  FGGE 
best  meets  this  requirement.  However,  the  study  should  also  use  concurrent  ocean 
observations  for  verification  of  the  mixed  layer  predictions.  Once  the  causes  of  the 
Arabian  Sea  mixed  layer  cooling  are  known,  then  the  relationship  of  this  cooling  to 
other  monsoon  characteristics,  such  as  rainfall,  may  be  more  precisely  analyzed. 
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APPENDIX  A 
THE  NAVY  OPERATIONAL  GLOBAL  ATMOSPHERIC  PREDICTION  SYSTEM 

The  Navy  Operational  Global  Atmospheric  Prediction  System  (NOGAPS)  at  the 
Fleet  Numerical  Oceanography  Center  (FNOC)  is  a  slightly  modified  version  of  the 
UCLA  general  circulation  model.  NOGAPS  has  been  the  Navy's  operational  atmos- 
pheric forecast  model  since  August  1982.  The  following  sections  describe  the  various 
features  of  NOGAPS  as  used  during  the  experiment.  This  information  is  adapted  from 
Ranelli  (1984).   The  complete  model  has  been  described  by  Rosmond  (1981). 

1.        DYNAMICS 

The  dynamics  of  the  UCLA  GCM  are  described  in  detail  by  Arakawa  and  Lamb 
(1977)  and  are  only  discussed  briefly  here.  NOGAPS  is  a  primitive  equation  model. 
The  prognostic  variables  are  horizontal  velocity,  V,  temperature,  T,  surface  pressure, 
ps,  and  specific  humidity,  q.  Additional  prognostic  variables  associated  with  the  plari- 
etary  boundary  layer  (PBL)  will  be  described  below.  The  finite  difference  scheme  has  a 
spatial  resolution  of  2.4°  lat.  by  3.0°  long.  The  variables  are  staggered  in  the  hori- 
zontal according  to  the  Arakawa  scheme  C.  The  center  grid  point  contains  the  T 
value.  The  meridional  wind  component,  v,  is  carried  at  points  north  and  south  of  the 
center  point  and  the  zonal  wind  component,  u,  is  carried  at  points  east  and  west  of  the 
center  point.  The  numerical  differencing  scheme  is  both  energy  and  enstrophy 
conserving. 

NOGAPS  uses  a  sigma  coordinate  system  in  the  vertical  defined  as: 

<T  =  (p  -  Pj)/7t 

where: 

p^  =  50  mb  and  7t  =  ps  -  p-  , 

p  is  pressure  and  p„  is  surface  pressure.  There  are  six  model  layers  in  the  vertical  with 
the  top  of  the  model  atmosphere  at  50  mb.  All  prognostic  variables  except  vertical 
velocity,  <y,  are  carried  at  the  middle  of  each  layer.  Vertical  velocity  is  carried  at  the 
layer  interfaces. 
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NOGAPS  uses  a  second  order  (leapfrog)  time  differencing  scheme  with  a  four 
minute  time  step.  Model  diabatics  are  executed  every  forty  minutes.  A  Matsuno  time 
step  is  used  every  fifth  time  step.  This  is  used  to  control  the  computational  mode  and 
to  assist  in  the  assimilation  of  the  diabatic  effects. 

2.        MODEL  DIABATICS 

The  sophisticated  model  diabatics  contained  in  NOGAPS  is  an  important 
component  in  this  experiment.  This  treatment  of  the  diabatic  processes  is  necessary  to 
adequately  simulate  fluxes  across  the  air-sea  interface  and  to  propagate  the  full  effect 
of  these  changes  throughout  the  atmosphere.  NOGAPS  directly  computes  the  physical 
processes  for: 

dry  convective  adjustment 

large-scale  precipitation 

diagnosis  of  stratus  cloud  depth 

mid-level  convection 

ground  hydrology 

surface  friction 

horizontal  diffusion  of  momentum 

radiative  transfer  processes 

cumulus  convection 
1.    Planetary  Boundary  Layer 

The  planetary  boundary  layer  (PBL)  is  defined  as  a  well  mixed  layer  in  mois- 
ture, moist  static  energy  and  momentum.  It  is  assumed  to  be  capped  by  discontinuities 
in  temperature,  moisture  and  momentum.  The  PBL  treatment  in  this  model  follows 
DeardorfT(1972)  and  has  been  formulated  for  the  UCLA  GCM  by  Randall  (1976).  It 
allows  for  interaction  between  the  PBL  and  cumulus  cloud  ensembles  and/ or  a  stratus 
cloud  layer  at  each  grid  point.  Surface  fluxes  are  determined  using  a  bulk  Richardson 
number  based  on  the  values  of  the  sea  surface  temperature  and  the  values  of  V,  T  and 
q  from  the  adiabatic  portion  of  the  model.  These  values  are  then  used  to  predict  a  new 
PBL  depth  and  the  strength  of  the  inversion  jumps. 

The  NOGAPS  PBL  is  constrained  to  remain  in  the  bottom  sigma  layer  of  the 
model.  This  differs  from  the  original  formulation  of  the  UCLA  GCM,  in  which  the 
PBL  was  allowed  to  pass  out  of  this  layer.  An  overly  deep  PBL  can  result  in  serious 
computational  problems  with  the  model.  Constraining  the  PBL  in  this  way  imposes  a 
maximum  depth  of  about  200  mb  on  the  PBL. 
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2.  Cumulus  Parameterization 

Cumulus  parameterization  in  NOGAPS  follows  the  scheme  of  Arakawa  and 
Schubert  (1974)  as  introduced  into  the  model  by  Lord  (1978).  In  the  model,  cumulus 
clouds  must  have  their  bases  at  the  top  of  the  PBL.  Cloud  tops  can  be  at  all  sigma- 
levels  above  the  PBL.  Cumulus  clouds  are  modeled  as  entraining  plumes  in  which 
environmental  air  is  mixed  with  the  PBL  air  from  which  the  cloud  originated. 
Tendencies  of  moisture,  temperature  and  momentum  are  diagnosed  as  well  as  the  cloud 
mass  flux.  The  cloud  base  mass  flux  removes  mass  from  the  PBL,  which  decreases  the 
PBL  depth.  Condensation  occurs  at  each  grid  point  where  the  air  becomes  supersatu- 
rated. A  moist  convective  adjustment  procedure  removes  convective  instability 
between  mid-tropospheric  layers  that  is  not  eliminated  by  clouds  originating  from  the 
PBL. 

3.  Radiation 

The  radiation  parameterization  follows  Katayama  (1974)  and  Schlesinger 
(1976).  It  includes  both  a  diurnal  variation  and  interaction  with  the  cloud  distribution. 
Radiative  transfer  processes  for  incoming  solar  radiation  are  computed.  Effects  of 
water  vapor,  Rayleigh  scattering  by  air  molecules  and  absorption  and  scattering  by 
water  droplets  in  clouds  are  included.  Reflection  due  to  clouds  is  also  calculated.  The 
model  cloud  cover  predicted  by  the  PBL,  the  cumulus  parameterization  and  large-scale 
precipitation  interact  with  the  long-wave  radiation.  The  net  surface  heat  flux  is 
computed  as  a  function  of  the  incoming  solar  heat  flux,  long-wave  radiation  flux, 
evaporative  heat  flux  and  sensible  heat  flux. 

3.        PROBLEMS  WITH  NOGAPS  IN  THE  TROPICS 

The  solar  radiation  flux  is  computed  indirectly  from  the  cloud  parameterization 
calculations  and  the  water  vapor  distribution.  The  solar  radiation  flux  analyses  in 
Chapter  II  indicate  that  this  scheme  is  not  appropriate  in  the  tropical  Arabian  Sea 
because  the  solar  radiation  flux  does  not  decrease  at  all  after  onset.  This  is  not  consis- 
tent with  the  findings  of  Gautier  (1986).  The  NOGAPS  parameterization  scheme  is 
more  appropriate  for  mid-latitude  regions.  Therefore,  NOGAPS  should  be  used  selec- 
tively in  tropical  regions. 
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APPENDIX  B 
DATA  HANDLING  PROCEDURES 

1.  ATMOSPHERIC  FORCING  FIELDS 

The  NOGAPS  and  NASA  atmospheric  forcing  fields  are  available  at  6-h  inter- 
vals. The  analyses  are  at  00  and  12  GMT  and  the  6-h  predictions  are  at  06  and  18 
GMT.  The  NOGAPS  atmospheric  forcing  fields  for  this  study  cover  the  area  from  the 
equator  to  25°N  and  from  50°E  to  75°E,  with  2.5°  grid  spacing.  The  NASA  atmos- 
pheric forcing  fields  cover  the  area  from  2°S  to  26°N  in  4°  intervals  and  from  50°E  to 
75°E  in  5°  intervals. 

The  NOGAPS  analyses  are  missing  data  at  all  grid  points  on  several  days  during 
the  onset  period.  The  missing  data  are  replaced  with  the  previous  24-h  values  of  that 
forcing  field  at  each  grid  point.  This  sets  the  missing  value  close  to  the  recent  trend 
value  instead  of  setting  it  to  an  average  value  that  might  be  far  different  from  the 
recent  trend. 

The  data  are  analyzed  in  three  ways.  First,  a  10-day  average  of  each  atmospheric 
forcing  field  at  each  grid  point  is  computed  for  the  pre-onset  period  ending  3-4  days 
before  onset.  Second,  ten-day  averages  are  computed  for  each  atmospheric  forcing 
field  for  the  post-onset  period  beginning  3-4  days  after  onset.  The  pre-onset  and  post- 
onset  ten-day  averages  are  computed  for  00,  06,  12  and  18  GMT,  although  only  the  12 
GMT  10-day  averages  are  reported  in  Chapter  II.  Finally,  a  time  series  of  each 
atmospheric  forcing  field  is  analyzed  by  using  a  5-grid  point  spatial  average  around 
each  of  the  four  analysis  points  in  Fig.  2.1.  The  time  series  begins  at  00  GMT  approx- 
imately eight  days  prior  to  onset  and  continues  until  18  GMT  approximately  ten  days 
after  the  first  day  of  onset. 

2.  MIXED  LAYER  MODEL  DATA 

The  Garwood  mixed  layer  model  requires  atmospheric  forcing  input  on  an  hourly 
basis.  This  means  that  the  NOGAPS  and  NASA  atmospheric  forcing  fields  must  be 
interpolated  since  they  are  analyzed  every  six  hours.  As  stated  in  Chapter  III,  a 
sophisticated  scheme  (Gallacher,  1979)  is  used  to  interpolate  the  solar  radiation  flux 
fields.  The  scheme  uses  the  location  of  the  grid  point  and  time  of  year  to  estimate  the 
diurnal  changes  in  solar  radiation  flux.  The  wind  stress  and  net  heat  flux  are  linearly 
interpolated  between  the  6-h  intervals. 
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